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Coexistence of antiferromagnetic order with superconductivity in many families of newly dis¬ 
covered iron-based superconductors has renewed interest to this old problem. Due to competition 
between the two types of order, one can expect appearance of the antiferromagnetism inside the 
cores of the vortices generated by the external magnetic field. The structure of a vortex in type 
II superconductors holds significant importance from the theoretical and the application points of 
view. Here we consider the internal vortex structure in a two-band s± superconductor near a spin- 
density-wave instability. We treat the problem in a completely self-consistent manner within the 
quasiclassical Eilenberger formalism. We study the structure of the s± superconducting order and 
magnetic field-induced spin-density-wave order near an isolated vortex. We examine the effect of 
this spin-density-wave state inside the vortex cores on the local density of states. 

PACS numbers: 74.20.De, 74.25.0p, 74.25.Ha 


I. INTRODUCTION 

The emergence of superconductivity at the onset of 
magnetism is a hallmark of many families of unconven¬ 
tional superconductors. Recently discovered iron based 
superconductors (FeSCs) provided a new addition to this 
list. Parent compounds for many of the FeSCs have 
the spin-density-wave (SDW) order, and superconduc¬ 
tivity (SC) appears upon doping or under pressure.^^lt^ 
The doping-temperature phase diagrams describing lo¬ 
cation of these two phases vary from material to ma¬ 
terial. In some systems, for example in the 1111 fam¬ 
ily RFeAsOi-jjFj, (R is a rare-earth element), the SDW 
phase abruptly disappears once the SC phase devel¬ 
ops. In other systems, such as 122 compounds based 
on BaFe 2 As 2 , the SDW order coexists with SC within 
some range of parameters. 

To understand how the SDW and SC phases interact 
with each other and what triggers superconductivity with 
such high transition temperatures, knowledge of the pre¬ 
cise structure of order parameters is essential. The local 
electronic structure near defects can be extracted from 
scanning tunneling spectroscopy (STM) measurements 
and it provides vital information about the order parame¬ 
ter. A Defect could be either an impurity or a topological 
singularity like a vortex induced by magnetic field. Here 
we focus on structure of an isolated vortex. 

It is well known that the shape of the vortex and 
the electronic structure close to it are very sensi¬ 
tive to the gap structure.!^ Copper oxide based high- 
temperature superconductors have been subjected to ex¬ 
tensive research for various kinds of competing orders 
inside the vortex cores. The signatures of such vortex- 
core orders have been reported in Bi 2 Sr 2 CaCu 208 -r<;P 
La2-5Sr5Cu04,l^ YBazCuaOT-i,^^ YBa2Cu408,l^ 
and Tl2Ba2Cu06-i-5p^ From the standpoint of the the¬ 
ory, several different approaches have been adopted to 
explain these experimental observations. Arovas et 
and Sachdev et studied antiferromagnetism in the 
vortex cores within phenomenological Ginzburg-Landau 


free-energy functional method. Ghosal et used mi¬ 
croscopic Bogoliubov-de Gennes (BdG) technique to in¬ 
vestigate this problem for the superconductors with a 
d-wave symmetric order parameter. The BdG technique 
was used heavily by many researchers to understand var¬ 
ious aspects of the competing orders inside the vortex 
cores.l23t^We also mention the work of Garkusha et al^ 
where the Usadel equation formalism has been used to 
explore the problem of antiferromagnetic vortex cores. 
The Usadel equations, however, are only applicable in 
the dirty limit, when the electronic mean free path is 
shorter than the coherence length, hence only appropri¬ 
ate for s-wave superconductors. 

The vortex state in FeSCs has been studied extensively 
by the STIVP^ and several novel features near vortex 
cores have been revealed. In the first study of the vor¬ 
tex structure in the optimally-doped BaFei, 8 Coo. 2 As 2 by 
Yin et aZ.^ no subgap states have been found. This is 
most probably due to large quasiparticle scattering rate 
in this material. On the other hand, optimally-doped 
Bao. 6 Ko. 4 Fe 2 As 2 does show peak at the vortex center, 
which is shifted from the Fermi level to lower energy.l^ 
This shift was attributed to the quantum effect, which 
is realized in the materials with moderate values of the 
product of the Fermi momentum kp and the coherence 
length ^0 at temperatures lower than Td(kp^o). Alter¬ 
natively, such energy shift of the localized state which 
breaks the particle-hole symmetry can be caused by mag¬ 
netic field-induced order in the vortex core. This sce¬ 
nario is very likely when a superconductor is close to a 
SDW instability. This possibility, however, has not been 
considered in Ref. 1341 Similar downshift was found in 
LiFeAs by Hanaguri et aZ.,1^ even though this material 
does not have obvious proximity to magnetism. Song 
et aZ.I^ studied the vortex state in FeSe and found en¬ 
hanced C 4 symmetry breaking in the vortex core which 
is probably related to orbital order in this material. 

These compelling features have motivated many the¬ 
oretical works. One class of theories has associated 
the particle-hole asymmetric finite-energy peaks to the 
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normal-state band structure of the materialsP^^^In this 
case the mechanism of particle-hole symmetry breaking is 
due to the quantum effect discussed in Ref. [351 Contrary 
to this proposal, several other authors have considered 
orbital^, nematid^ or SDW ordeil^^HllI, Hung et 
have included orbital ordering within a self-consistent 
BdG approach, and explained the enhanced C 4 symme¬ 
try breaking observed in FeSe by Song et al. in Ref. 1331 
Similar results were reported by Jiang et and Hu et 
al^ for the SDW order also using the BdG method. 

In this paper, we consider the emergence of the SDW 
order in the vortex cores and its spectroscopic conse¬ 
quences. We use the quasiclassical Eilenberger approach 
to study the problem of the field-induced SDW order in¬ 
side an isolated vortex. Both the BdG and Eilenberger 
approaches have their own advantages and complement 
each other. The BdG method is more microscopic. On 
the other hand, the Eilenberger approach relies on few 
most essential physical parameters. It is numerically 
less expensive and allows to study more complex prob¬ 
lems. We compute distribution of the superconducting 
and SDW order parameters inside the core and typical 
length scales for both order parameters. We also investi¬ 
gate influence of emerging SDW order on the density of 
states (DOS) near the vortex. This paper is organized in 
the following manner. In the next Sec. |TT] we describe 
the details of the model and the method. In Sec. m we 
discuss the results and conclude in Sec. IIVI 


II. MODEL & METHOD 


A. Quasiclassical equations for a two-band 
superconductor with spin-density wave 


We consider a simple minimal model with two cylindri¬ 
cal Fermi surfaces, which allows us to capture qualitative 
understanding of the problem. For the dispersion of the 
holelike Fermi surface, we take 


a(k)^a(k)=/..-;^, (1) 

2mh 

and for the electronlike Fermi surface we consider follow¬ 
ing dispersion. 


Ce(k) = 6(k-Q) 


(kx-Qxf 

2 TOe(l-e) 


{ky Qy) 

2 me (1-|-e) 


Me: (2) 


where {Qx,Qy) is the SDW ordering vector and /r/i:Me 
are the energy offset for the hole and the electron band 
respectively. Fig. shows a schematic picture of the two 
Fermi surfaces. It is useful to write these dispersions as. 


4(k) = -C, (3) 

U^ + Q)=^ + 2S, (4) 

where <5 is the energy scale, which measures the deviation 
from perfect nesting. In general, d is a function of the 
angle on the Fermi surface (j) and goes to zero at the hot 



FIG. 1. (Color online) Schematic representation of the hole- 
like and electronlike Fermi surfaces (solid lines) centered 
around the F point and M point respectively. A shifted hole- 
like Fermi surface is shown with dashed line. Filled circles are 
the hot spots, where the nesting is perfect. 


spots (shown in Fig[^. For the dispersions considered 
here, 

“ ^iso T *^ani COS 20, (5) 


with 


^iso — .- 
^ani — 


rrihUh 


2 Vme(l - e^) 
ruhlJ-h e 


2me 1 — 

and we treat Siso and dani as tuning parameters. 

The model Hamiltonian is same as used by several 
other group^^^®^. 


H — Hkin T Use T HgUw; 
Hkin = 'y ' ^i(k)c| k,aG,k,ctj 


Hsc = ^ A, (*u^)„^ c) k^„cI.-k./3 + 


2,k,Q; 

Ja/3 

2,k,Q;,/3 

Hsdw = (a cj ^^C 2 ,k ,,9 + h.c. 

k,a,/9 


( 6 ) 

( 7 ) 

, ( 8 ) 
( 9 ) 


where are the SG order parameters for two bands 
with f = 1,2 being the band index and M is the SDW 
order parameter. We only consider singlet superconduc¬ 
tivity. For incommensurate SDW order M is a complex 
quantity. Here we consider only the commensurate SDW 
order, which makes M a real quantity. We will briefly 
discuss the consequences of incommensurability in the 
SDW order. The indices a, /3 denote the spin states, 
and c]^ jj. ^ (cqk,a) is the fermionic creation (annihilation) 
operator for a fermion in the band with spin a. We 
consider s± state for superconductivity with equal gap 
magnitudes in two bands with a relative sign change. 

The self-consistency conditions read, 

(Cj -k,aCj.k./3) , (10) 


( 11 ) 
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Here and are the pairing interactions for the 

SC and SDW phases respectively and assumed to be mo¬ 
mentum independent. 

In the extended particle-hole basis, = 

(cj ij. ci.-k,^, C 2 C 2 ,-k,^), the Hamiltonian reads, 


H = ^ 

k 

■ Cl Ai M 0 ■ 

AJ -Cl 0 M 

M 0 C 2 A 2 

.0 M -C 2 . 


( 12 ) 

(13) 


The 4x4 matrix Green’s function for this mean-field 
Hamiltonian is 


g = (ioji -n) \ (14) 

where ui = 2TTT{n+ 1/2) is the fermionic Matsubara fre¬ 
quency and 1 is 4x4 identity matrix in the two band 
particle-hole space. 

Next, we derive the quasiclassical equations, which 
were first obt ained by Eilenberger for conventional 
superconductor^^HEll^ These are transportlike equations 
for the kinetic energy integrated Green’s functions, 


9 = 


i 

TT 


rfC7 • 


(15) 


where 7 is a 4x4 diagonal matrix with elements 
(1,-1,—1,1). In compact matrix form the quasiclassi¬ 
cal equation reads. 


[uj + —vf • a) 7,g 
\ ic / . 


vf • Vg 


{Hs + iisc + 'Hsdw)7 ,9 


= 0 , 


(16) 


where vp is the Fermi velocity and A is the vector po¬ 
tential. Hsc, 'Hsdvj are the SG and the SDW components 
of the mean-held Hamiltonian in the basis spanned by di. 
TLs is a 4x4 diagonal matrix with elements (0,0,25,-2(5) 
containing information about nesting between the Fermi 
surfaces. Its contribution drops out in the absence of 
the SDW order. This equation agrees with one derived 
by Moor et In the pure superconducting limit Eq. 

(16 1 reduces to the well-known Eilenberger equation. 


VF 


Vg-I- ^1 


■ Aj 7,5 


'Hsc1,9 


= 0. (17) 


Equation (16 1 has to be supplemented with the normal¬ 


ization condition. 


9^ — 9 bulk- (18) 

In particular, g^^ip. = 1 for a uniform superconductor 
without SDW order. The self-consistency conditions for 



FIG. 2. (Color online) Coordinate systems used to solve the 
Eilenberger equations. The real space lab frame is shown with 
solid lines, while the dashed lines represent the new coordi¬ 
nate system. Filled circle is the location of the vortex core 
and chosen as the origin. Open circle is the point where the 
solution is required. A point in the lab frame r = XpX -f yp-y 
maps to \pV + HpU in the new coordinate system. 


the order parameters can be expressed in terms of the 
Eilenberger functions as. 


Ai = ^ (E-gi2 + V,^^g34)F.s. - (19) 

UJ 

A2 = ittT {Vi2 9i2 + 1^22 534)f.s. - (20) 

UJ 

_^ "^sciw 

M = iirT E- (513 + 524 ) F.S. ■ ( 21 ) 


Here gij are components of 4 x 4 matrix Green’s function 
(indices i,j = (1, 2)/(3,4) correspond to the hole/electron 
band), {gij)p s angular average over the Fermi 

surface for the respective bands weighted with the den¬ 
sity of states, which is approximately the same for the 
both bands for the Fermi surfaces considered here. As 
for the pairing-interaction matrix, to get the s± super¬ 
conducting state we take into account only the interband 
repulsive interaction and neglect intraband terms, i.e., 
I’Ti = ^22 — 0- For pure SG state, Eq. (17) equations 
can be transformed to a set of Riccati equations, which 
makes the numerical solution much easier .^However this 
transformation is not useful for the problem considered 
here. In the following section, we discuss the strategy to 
solve these equations for the present case. 
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FIG. 3. (Color online) The phase diagram in the (T,5ani) 
plane for the two phases at zero magnetic field. Thick 
black/red line indicates transition to the SDW/SC state. The 
thin dashed line shows the Ts in the absence of the SC corre¬ 
lations. Filled squares are considered as representative cases 
in this paper. 


B. Numerical solution 

The Eilenberger equations are first-order partial differ¬ 
ential equations. A standard tool for solution of this type 
of equations is the method of characteristics. The basic 
idea of this method is to introduce the new coordinate 
system, in which the partial differential equation reduces 
to an ordinary differential equation. Hence, it is useful 
to introduce a coordinate system spanned by the two or¬ 
thogonal vectors, a unit vector along the direction of the 
Fermi velocity v and a unit vector u orthogonal to v, see 
Fig. § The unit vectors spanning the new coordinate 
system read, 

D = cos0i-I-sin0y, (22) 

u =—sinOx-\-cosOy. (23) 

Here 9 is the angle between the Fermi velocity and the x 
axis in the lab frame. A point in the lab frame r = (cc, y) 
maps to (A, ^) in the v-u frame. A point (xp, yp) at which 
solution is desired transforms to {Xp, yip) as, 

Xp = Xp cos 9 + yp sin 9, (24) 

fip = —Xp sin 9 + yp cos 9, (25) 

where the parameter fip has the meaning of an impact 
parameter. For a fixed trajectory, this impact parameter 
is uniquely determined by {Xp, yp) and it does not change 
with change of A. The quasiclassical equations are solved 
along these classical trajectories (u) in the real space. 
Along such trajectories quasiclassical equations reduce to 
system of ordinary differential equations, which are much 
easier to handle than solving a set of partial differential 
equations. Far away from the bulk, the system is ho¬ 
mogeneous. The homogeneous values are used as initial 


values. Now on a given trajectory, there are two possi¬ 
bilities. One can integrate towards the defect (vortex in 
this case) from the two extreme ends (A = ±oo) of the 
trajectory. Due to the first-order nature of the equations, 
the numerical solution readily converges to exponentially 
growing function. Of course, these exponentially growing 
solutions are unphysical. However, it is possible to con¬ 
struct the physically bounded solution at any point using 
the exponentially growing solution using the explosion 
method. The explosion method exploits exponentially 
growing solutions to obtain the physical solution.[S3t^ 
(See Appendix for details) For each point in the real 
space, one has to solve the Eilenberger equations for all 
the trajectories and for each Matsubara frequency. To 
obtain a physical solution, we solve the Eilenberger equa¬ 
tions from two opposite directions A = ±oo towards the 
point, where the solution is desired. As shown in the Ap¬ 
pendix, the two exploding solutions g± diverging in the 
±oo limits, provide the physical solution gp as 


9-9+ - 9+g- 
9-9+ + 9+9-' 


(26) 


Once all the Eilenberger functions are computed for an 
initial guess for the order parameters, an updated set of 
order parameters is recalculated, and this process con¬ 
tinues till it converges to a solution. It is convenient to 
normalize all the energy scales to and all the lengths 
are measured in the unit of superconducting coherence 
length ^0 = vf/{2ttTc). Here vp '^s the average Fermi ve¬ 
locity of the two bands. We consider weak ellipticity for 
the electronlike Fermi surface, and the Fermi velocities 
of the two bands are roughly equal. All our results are 
presented in these units. 


III. RESULTS & DISCUSSION 

Coexistence of the SDW state and the superconduc¬ 
tivity is very sensitive to the underlying electronic struc¬ 
ture. For the two-band model we consider, the nesting 
function (5(0) in Eq. (§ can be tuned to get a co-existing 
phase.SSHm Here we are interested in a situation where 
there is no long-range SDW order in the absence of the 
magnetic field. Fig. |^shows the phase diagram as a func¬ 
tion of the anisotropic nesting parameter (5ani for a hxed 
value of (5iso = 0.16(27rrso)- The presence of the super¬ 
conductivity strongly modihes the SDW phase bound¬ 
ary (thin dashed line in Fig. [^. The region between 
the original and SC-renormalized phase boundaries pro¬ 
vides a possibility of the SDW order in the regions, where 
the SC phase get suppressed locally. The phase diagram 
shown in Fig. [^is only includes the commensurate SDW 
phase. Vorontsov et have shown that the incom¬ 
mensurate SDW phase may co-exist with the SC in a 
larger area of the phase diagram. We consider few rep¬ 
resentative cases with (5ani/27rTso = 0.25,0.26 and 0.27. 
Here T^q is the SDW transition temperature for a sys¬ 
tem with perfect nesting and with the same interaction 
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FIG. 4. (Color online) The temperatnre dependence of the 
bulk superconducting order parameter and the SDW orders 
at the center of the vortex for three values of the parame¬ 
ter 5ani/27rrso= 0.25,0.26 and 0.27. All the energy scales are 
normalized to T^. The thin line shows the temperature de¬ 
pendence of the magnitude of the SC order in the bulk, which 
is almost identical for two bands. This signs of SC order pa¬ 
rameters are opposite for two bands. The mean-field SDW 
transition temperature in the absence of superconductivity is 
indicated with filled square for each case. 



FIG. 5. (Color online) The spatial variation of the magni¬ 
tude of the SDW order at T = 0.3Tc for 5ani = 0.25(27rTBo), 
where is the SDW transition temperature for a system 
with perfect nesting and with the same interaction strength. 
The magnitude of the SDW order is normalized to Tc. The 
hot spots are located near the y-axis. 


strength. We have set T^o = 2Tco, where Tco is the su¬ 
perconducting temperature. As we are mostly interested 
in low-temperature behavior, we restrict our calculations 
below T = O.STj,. 


A. Order parameter profiles 

Fig.i shows the temperature dependence of the bulk 
superconducting gap and the SDW order in the core. It 
is evident that the SDW order appears roughly below 
the temperature one would expect from the phase di¬ 
agram shown in Fig. The temperature dependence 
of the SDW order parame ter deviat es strongly from the 
mean-field behavior (oc a/1 — T/Tc). We see two differ¬ 
ent temperature regimes. At lower temperature the SDW 
order grows strongly, but slightly below the mean-field 
SDW transition Tg it develops a tail, which survives even 
above Tg. It should be noted here that the phase diagram 
is based on the commensurate SDW phase for the normal 
state electronic structure. However, the electronic struc¬ 
ture of the vortex core states is not the same as in the 
normal state. The onset of the SDW order is mostly de¬ 
termined by the core bound states. The SDW transition 
temperature T^ cx exp[-l/lVoFsdw/(<5iso, i^ani)], where 
the function /((5iso, i^ani) depends on nesting parameters, 
I4dw is the SDW interaction and Nq is the density of 
states. As in the vortex core the density of states is 
higher than the normal-state value, the SDW onset tem¬ 
perature may exceed the mean-field value. 

We restrict ourselves to the commensurate case only, 
but for the incommensurate case the phase boundary 
shifts towards slightly higher temperatures. As reported 
by Vorontsov et al. in Ref. 1461 an incommensurate or¬ 
der may exist in a larger portion of the phase diagram. 
We have also performed calculations, where we allow in¬ 
commensurability in the SDW order. With incommen¬ 
surability the SDW order parameter becomes complex 
and acquires a finite phase. We found that the phase 
of the SDW order parameter is temperature dependent 
and varies very weakly in the real space. The incommen¬ 
surate order persists above the phase boundary shown 
in Fig. 1^ Since we did not find anything qualitatively 
different, we will focus on the commensurate case only. 
Furthermore, there is no qualitative difference between 
the cases considered here, except for the temperature de¬ 
pendence. Therefore, we continue our discussion with 
(5ani/27rTso = 0.25. For this Jani and (5iso/27rTso = 0.16 
the nesting hot-spot angle (f) « 64.9° is close to the y di¬ 
rection which strongly influences anisotropic properties 
of the vortex. Figure shows the magnetic field-induced 
SDW order parameter in the real space at T = 0.3Tc. 
Spatial coordinates have been normalized to the super¬ 
conducting coherence length and the SDW order is 
normalized to T^. An important feature is the oscilla¬ 
tions of the SDW order along the x direction which is 
most clearly seen in the inset. At lower temperatures 
when the SC vortex is very small, the SDW state is lo¬ 
calized very close to the vortex core. As the tempera¬ 
ture grows, vortex becomes larger and the region with 
the SDW order also increases, due to a larger region of 
the suppressed superconductivity. The size of the vortex 
is larger in the presence of the field-induced SDW order. 
This property can be seen in Fig. in which we compare 
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FIG. 6. (Color online) The left panel shows the gap magnitude normalized to its bulk value at T = O.STc in the presence of 
the SDW order and the right panel shows the gap structure at the same temperature with no SDW order for the hole band. 


a SC vortex with and without the magnetic field-induced 
SDW order. The SDW order makes vortex larger and 
anisotropic. The intrinsic anisotropy of the underlying 
band structure is weak. Hence the large anisotropy in 
the real space is mainly due to the field-induced SDW 
order. Strong enhancement of the anisotropy is reflected 
in the characteristic length scales associated with the SC 
order along the two principal directions. Fig. [^shows the 
length scales associated with the SDW and the SC 

order along the x and y axis. We define the super¬ 

conducting coherence length as a distance from the 
vortex core, where the order parameters reaches half of 
its bulk value. Similarly, the magnitude of the SDW or¬ 
der drops to half of its value at the core at a distance 
^sdw from the vortex core. As illustrated in Fig. the 
SC length scales along x and y directions become differ¬ 
ent in the presence of the SDW order and this is mainly 
due to anisotropy in the field-induced SDW order. Note 
that the SDW correlations are stronger along the y di¬ 
rection which is closer to the nesting hot spots. This 
causes stronger suppression of the SC order and reduces 
the SC characteristic length along this direction. Next, 
we discuss the density of states near the vortex core. 


B. Density of states 

All the previous calculations were done using the Mat- 
subara frequencies. For the DOS calculation it is nec¬ 
essary to go to the real frequencies. We used the order 
parameter profiles calculated in the Matsubara represen¬ 
tation. The analytic continuation, iuj ^ E + iy, is done 
with an artificial broadening rj = 0.057);. This process 
also requires solution of the Eilenberger equations for 
the real frequencies using the same approach as in the 
calculation of the order parameters in the previous sec¬ 
tion. The total DOS is the sum of the partial DOSs for 



FIG. 7. (Color online) The temperature dependence of char¬ 
acteristic length scale for SDW order (^Jy^) and the charac¬ 
teristic length scale of the SC order (^J/j,) on the hole band 
for (Jani/SyrTso = 0.25. Subscripts denotes x and y directions 
in the real space. The same characteristic length scale for a 
pure superconductor is plotted with a dotted dashed line for 
comparison. 

each band and, in terms of the Eilenberger functions, it 
is given by 

N{E) = Re{[gii{E + iy) + gs^iE + ig)])FS- (27) 

Fig. i presents evolutions of the density of states along 
the two principal directions for two temperatures. The 
first row of Eig. shows DOS along the x axis, which is 
away from the hot spots and the second row shows DOS 
along the y axis which is closer to the hot spots. The 
field-induced SDW order enhances the violation of the 
C 4 rotational symmetry the vortex center. Eor conven¬ 
tional superconductors, the DOS is always particle-hole 
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FIG. 8 . (Color online) DOS along the x direction (first row) and along y direction (second row) at T^O.STc in panel (a), (b) 
and at T=0.5Tc in panel (c) and (d) for 5ani/‘^'xTc=Q.2&. 


symmetric, unless the superconductor is in the quantum 
regime when is not too large and T <C Tc/kp^o^ 
Another key feature of the classical clean-limit DOS in 
the vortex core is the sharp peak at zero energy corre¬ 
sponding to the localized state. The quantum effects shift 
this zero-bias peak to the finite energy but, do not break 
the rotational symmetry. The emergence of the SDW 
order leads to particle-hole asymmetry in the DOS and 
also strongly violates the C 4 symmetry. The particle-hole 
asymmetry is bigger along x direction because of stronger 
deviation from nesting in this direction. Another impor¬ 
tant feature which is visible in the Fig. is suppressed 
spectral weight at the core indicating that the energy of 
localized state is shifted from zero to a finite value cor¬ 
responding to opening of a minigap in the core. This is 
shown more clearly in Fig. [^in which we show the DOS 
plots at several representative points. This small gap in 
the DOS near the vortex core is particle-hole asymmetric. 


which is a hallmark of energy gap due to the SDW order. 
This gap vanishes away from the vortex core indicating 
presence of a state with energy close to zero localized 
outside the core. As the temperature increases and the 
SDW order weakens, the apparent gap in the core dis¬ 
appears as shown in panel (c) and (d) of Fig. |^at T = 
0.5Tc. The described features are the keys to distinguish 
between the quantum effect and the field-induced SDW 
order. Figure also shows the DOS for two different 
temperatures. The temperature dependence of the DOS 
is easy to understand. As the temperature increases, the 
SDW order weakens, which reduces the degree of C 4 sym¬ 
metry breaking and the particle-hole asymmetry in the 
DOS. 





































FIG. 9. (Color online) Panels (a) and (b) show DOS near the 

IV. SUMMARY AND CONCLUSION 


We study the structure of an isolated superconducting 
vortex near a SDW instability inside the superconduct¬ 
ing dome. We show that the SDW order develops inside 
the vortex below a critical temperature determined by 
the strength of the SDW instability. This leads to C 4 
symmetry breaking near the vortex core. If there is al¬ 
ready C 4 breaking in underlying band structure, then it 
gets strongly enhanced due to the SDW order near the 
vortex core. The corresponding deformation of the vor¬ 
tex shape can be imaged by the STM technique. We 
find that the field-induced SDW order persists beyond 
the superconducting vortex region. The tunneling DOS 
carries very strong signatures of this field-induced order. 
A small energy gap develops inside the core and gives rise 
to strong particle hole asymmetry, which is pronounced 
along the directions away from the hot spots. Our re¬ 
sults are in qualitative agreement with the STM data on 
on Bao, 6 Ko. 4 Fe 2 As 2 p^ which may indicate the presence of 
the vortex-core SDW order in this material. Our findings 
also agree with BdG-based works by other groups. 
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vortex core along x and y directions respectively at T^O.STc. 

Appendix A: Explosion method 

The general structure of the Eilenberger equations is, 

iur.jl. (Ai) 

which is a first-order ordinary differential equation. It 

is straightforward to show that if 5 is a solution of this 
equation then if' is also a solution, which implies 

= ao + Oi 5 , (A2) 

where oq is a constant matrix and ai is a complex num¬ 

ber. It can be further shown that the product of any two 
solutions of the Eilenberger equations is also a solution 
of these equations. This gives a very powerful relation, 

5^ = Oo = Qhulk^ (A3) 

which is very useful in obtaining the numerical solu¬ 
tions of these equations. There are multiple solutions 
to this system of equations. In pure superconducting 
state, there are three independent solutions. There are 
two divergent solutions along with a bounded physical 
solution. All the unphysical solutions decay to zero in 
the bulk. Let’s consider two such unphysical solutions, 
g± (X It can be shown that commutator of these 

two unphysical solutions gives the physically bounded so¬ 
lution, 


A± = g-g+±g+g-, 

(A4) 

A± = [r,A±]. 

(A5) 

The bounded physical solution is. 


II 

><> 

1 

(A 6 ) 
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the constant Cp can be determined using Eq. (A3) and 
(A 6 ) and it reads, 


1 



(A7) 


We use these unphysical solutions g+- in the bulk and 
integrate towards the vortex core starting from the bulk. 
Since these solutions grow exponentially, they can be 
easily computed numerically with appropriate boundary 
conditions. Far away from defects, we can ignore the spa¬ 
tial dependence of the order parameters. Since there is no 
long-range SDW order, we have the standard Eilenberger 
equations for the pure superconducting state in the bulk 
and in the basis we consider here, the Eilenberger Green’s 
function is a block diagonal matrix. The two bands are 
coupled through the self-consistency condition. There¬ 
fore it is sufficient to illustrate the idea for one band, for 
which we write down the equations explicitly, 

g = z{A*f + Af), (A 8 ) 

/ = - 2 iof - 2 iAg, (A9) 

P = 2ujp - 2iA*g. (AlO) 


We first find two unphysical solutions, which can be de¬ 
termined easily. Let’s consider, 

g = cie^^, (All) 

/ = C 2 e^^ (A12) 

P = cse^^. (A13) 

where = ±z/. Normalization condition requires, 

5 ± = 0. (A14) 


This ensures that all the unphysical solution decay to 
zero in the bulk. Which gives. 


cl + C 2 C 3 = 0, (A15) 

Cl = i^/c^c-i- (A16) 

Using these conditions, 


(C -I- 2w)c2 = 2 APC 2 C 3 , (A17) 

(C-2u;)c3 = 2A*Vh^h^. (A18) 

Thes e two e quations give the value of C = 
±. 2 p+ \A\'^/vf and 


C2 

C3 


-2iA 


Cl 


-2iA‘ 

(^vf — 2uj 


Cl. 


(A19) 

(A20) 


Here we fix ci = 1 and write the exploding solutions, 


g+ = exp 

g- = exp 

1 

P± = 


2QX' 

1 

-iAp+ 

+ 

Vf 

iA*p_ 

-1 

2QX 

1 

—iAp- 

'^F . 

. ^^*P+ 

-1 


ljj zL Q’’ 

Q = puj'^ + |A|2. 

Now we can write down the physical solution. 


9p = 


X_ 


X+ y/w 2 + |A |2 


uj —iA 
lA* —oj 


(A21) 

(A22) 

(A23) 

(A24) 

(A25) 


Once we get the values of these two diverging solutions in 
bulk then we can use the bulk values to integrate towards 
the vortex core and find the physical solution using two 
unphysical solutions. 
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